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Abstract 

Some puzzles which arise in matrix models with multiple cuts are pre- 
sented. They are present in the smoothed eigenvalue correlators of these 
models. First a method is described to calculate smoothed eigenvalue corre- 
lators in random matrix models with eigenvalues distributed in a single-cut, 
previous know results are reproduced. The method is extended to sym- 
metric two-cut random matrix models. The correlators are written in a 
form suitable for application to mesoscopic systems. Connections are made 
with the smooth correlators derived using the Orthogonal Polynomial (OP) 
method. A few interesting observations are made regarding even and odd 
density-density correlators and cross-over correlators in Z2 symmetric ran- 
dom matrix models. 



1 Introduction 



Matrix models have been used in a wide variety of applications, starting 
from quantum chaotic systems to condensed matter, QCD and string the- 
ory. The recent period has seen a large increase in our understanding of 
the properties of these models. In this work we have been interested in 
highlighting some unusual properties of two-cut random matrix models that 
have arisen in our study. The results are unexpected as they are not seen 
in matrix models when the density of eigenvalues has a connected support. 
Indeed there it is well-known [0, |2| that the correlator is universal, i.e. inde- 
pendent of the specific potential V which defines the probability measure. 
This is the basis for the theory under the universality of conductance fluc- 
tuations in mesoscopic systems ||3|. At first sight one is tempted to think 
that this universality persists when the potential is such that the support 
splits into disconnected segments. But it is found that, if indeed it is again 
universal, it belongs to a different universality class. If the standard large 
N-limit (the random matrices are x N) yields the smoothed correlation 
functions up here to an arbitrary constant, different methods report different 
results for this constant. Furthermore there are differences between these 
correlators when the size N of the matrices is an even or an odd integer. It 
is a rather intriguing phenomenon and, for instance, it is not clear how the 
naive renormalization group approach Q which consisted of integrating out 
one line and one row, could deal with such situations. We attempt here to 
understand and to give a unified picture of these results. 

The paper is divided as follows. It starts by establishing the notation 
and conventions and describes completely the method used for the model 
with a single-cut density of eigenvalues. Previously known results are repro- 
duced 1^, ^ . Then the method is extended to the model with two-cuts in 
the density of eigenvalues, restricted to symmetric potentials. Afterwards 
we devlop the formalism to include asymmetric potentials. Here an arbi- 
trariness remains as the constraint on the filling factor of the two parts of 
the support is not fixed at leading order in the large N-limit. The large 
N-equations for the correlator leave us with an undetermined constant C. 
Previous methods using the orthogonal polynomials and loop equations give 
different results for this constant . The orthogonal polynomial method is 
briefly outlined and the resulting correlators are sensitive to the even and 
oddness of the number of eigenvalues. Further, the constant C is different for 
the even and odd correlators found by the orthogonal polynomial method 
and that found by the loop equations. The conclusion summarizes these 
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results and attempts to give an explanation of these puzzles. 



2 Notation, Conventions 

We establish the notations and conventions and develop a method, which 
we extend to the two-cut model, to derive eigenvalue correlators for random 
matrix models with a single cut density of eigenvalues. 

Let us work with the unitary invariant ensemble of random Hermitian 
N X N matrices, with a probability distribution 

P{M) = ^exp{-NTrV{M)). (2.1) 



Define the operator for the density of eigenvalues 

1 ^ 

Pi^) = -jY.^{^-^^) (2-2) 

and 



i=i 



p{x) =< p{x) >= J P{M)^Tr6ix - M)d^^M. (2.3) 
Since P{M) gives a factor exp{—N'^ J V {x)p{x)dx) we have 

5p{x) 



5V{y) 

In the large N limit, we know that 



N' < p{x)p{y) >c . (2.4) 



' P^dy = \v'{x). (2.5) 
la X -y 2 

The solution is found through the averaged resolvent 

^ ^ N z-M Ja z-y ^ 

= W{z)-P{z)^{z-a){z-h). (2.6) 

Then 

^p(A) = P{X)^{X-a){b-\). (2.7) 

We then need to express P as a functional of V. There are various equivalent 
expressions, and in the following several will be needed. We follow here P] 
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and begin by multiplying eq. (p. 61) by ^"^jf — and integrate z over a 
large circle C, see Fig. 1. 

Since G{z) « ^ at infinity 

(i)- 



G{z)y^{z- a){z-b) dz _ ^ ^^^^ 

z — U 2i7T 



P{z){^{z-a){z-h)f dz 



z — U 2iTT 



P{u){u- a){u-h) (2.9) 



{in). 



V'{z)^{z- a){z-b) dz 
z — u 2m 



Therefore we obtain from eq. ( p.q ) 
1 = ]^V'{u)^{u-a){u-h) 



V'{u)J{u-a){u-h) 



1 , V'(x)Jix - a)ih- x) 
— / dx 



TT 



X — U 



(2.10) 



y'(x)V(a;-a)(6-x) 



X — u 



P{u){u- a){b-u). (2.11) 



Let u approach the real axis on the cut. The integral in eq. (2.11) has an 
imaginary part which cancels the first term of the right hand side which is 
pure imaginary. The real part of eq. ( |2.1lD gives 



P(A)(A-a)(6- A) = 1 + 



27r 



^ , V'{x)J(x-a)(b-x) 

dx 

X — X 



(2.12) 



I.e. 



P(A) 



vr V(A-a)(6-A) 



[1 + 



1 f\ V'{x)^{x-a){b-x)^ 



27r 



dx 



x — X 



(2.13) 
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Now we vary V, let us first ignore the variation of a and b (it is proved to 
be right below). Then 

^P(A) _ .dp{Xl ^.5p(A) Sa , ,9p(A), 6b 
<5y(/x) " ^5F(/x)^'^'' + ^ 9a '^''SVifi)^^ db ^^'''6V{ti)^ ' 

on the right hand side (r.h.s.) p is being treated as a function of V, a, b as 
given on the r.h.s. of eq. (|2.13|) . We show later that (|f )y,f) = 0. [{^)v,b is 
of course not the total derivative of p w.r.t. a]. Then 



6pW ^ 1 1 d ^{p-a){b-p) 

6V{pl) 27r2 ^(A-a)(6- A) dp X-p 

(2.15) 

and one verifies easily that the result is symmetric under exchange of A 
and p as it should. Note that the potential V has disappeared from the 
correlator, except indirectly through the end points a and b of the cut. This 
universality follows here trivially from the linearity of the {p, V) relation. 
Terms of the type 

have been ignored. The claim is that they vanish, but that's the only 
(slightly) tricky part. The representation eq. ( ^.13 ) is appropriate, among 



several other possibilities, because if one differentiates inside the integral 
with respect to a, it is still a meaningful integral. So let us calculate 

,dp, 1 1 1 1 



'da'''" 2(A-a)^' ' 2Tr^ ^(A-a)(6-A) 

If'' , Ib-x 1 



To prove that this is zero, let us return to eq. (|2.6|) and multiply it by 
^—^—^ and integrate again over a circle of large radius. Then 

z—a z—u o o o 



^,,,2; — 61, f dz 

G{z)\l- — -dz = <j) — = (2.18) 



z — a z — u Jc z 



as for large z, Giz) = \ ~ 1 and ~ -. While the second and 

° ' ^ ' 2 ' V Z—a z~u z 

third terms become 

.:^l£)(i_i:)ii = -P(„)(„-6) (2.19) 

c [z — u) Im 
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and 



Jc^ z — a z — u 2i7r 2\j u — a 2irJa^x — ax — u 

(2.20) 

Taking u = A + ze and using -^^^ = + ^7r5(a), the integral in eq. ( p.20| ) 
has a part which cancels the first term leaving 



2 Jc\ z — a z — u 2i7r 2tt Ja \ x — a x — u 

Repeating all the steps which led to eq. ( p.l5| ) i.e. combining eq. ( p.l8| ), eq. 
( 2.19| ), eq. ( 2.21| ) one finds an expression for p from 



^r,^^^. = P(A)(^-A) (2.22) 

2tt Ja \ X — a X — X 



which is 



p{\) 1 1 b-x 1 



/ r-,V'{x)dx (2.23) 

Ja y X — a (x — X) 



(A -a) 27rV(A - a){b -X)2Ja ^ x-a{x-X) 

' dp 
^ da 



thus proving that (f^)y,6 = 0. This completes the proof for the single-cut 



correlator. 



3 The Double Well 

Now let us extend the result to eigenvalues distributed in two disjoint bands 
([—a, —b] U [a, 6]). Let us first restrict ourselves to even potentials i.e. 

P(M) = Z'^exp{-NTrV{M)) 
P{-M) = P{M) (3.1) 

which implies for the resolvent 

G{-z) = -G{z). (3.2) 

Since we restrict ourselves to even potentials we cannot take a functional 
derivative of p{X) with respect to an arbitrary V{p), but we can fold the inte- 
grations over the positive part of the spectrum and then vary the potential. 
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Now 

TrV{M) = N I dXp{X)V{\) 



— oo 
oo 



/ dXV{XMX)+pi-X)]. (3.3) 







Consequently 

1 ^P(A) 
iV2 5V{p) 

where use has been made of 



< /'(A)p(m) >c + < p(A)/>(-/u) >c, (3.4) 



5y(/x) 

In the large N limit again 



^^'(")-5'(x-^). (3.5) 



2G{z) = V'{z)-P{z)^) (3.6) 

with (t{z) = (z^ — a2)(z2 — 6^). Note that this equation determines uniquely 
P{z), a k, h \ indeed take 

degV = 2n 

deg[P] = 2n - 3; 



P{z) has to be odd 

P{z) = aiz + a2Z^ + ■■■ + a„_iz^""^ (3.7) 

we thus have (n — 1) + 2 unknowns. Since G{z) ~2^oo - we have to fix 
the coefficients of eq. ( |3.6| ) at infinity from z^"""*^, z^""^, ■ ■ ■ , z^ , z~^ —>■ (n + 
1) conditions. Therefore no "filling" parameter creeps into the problem 
(although the question of spontaneous symmetry maybe eliminated by the 
assumptions here). 



Now we take eq. (^^), multiply by ^^"^2) integrate over a large 



circle in the z plane (using y ^^-u^ ^l^o been checked to give the same 
equation), see Fig. 2. We obtain 



2 ^ V\u)^^ 2i f-V'{x)./\^\ 



u 2iiT Jb x{x — u) 



^ 2i r~'' V'ix)./WixT\ P{u)a{u) 



2m J -a x{x — u) 



6 



which simphfies to 

P{u)a{u) V'{u)^/\a{u) 



2 + 



u 



u 



- -/ 

TT Jb 



1 



(x — u) (x + ii) 



We now let « = A + ie approach the cut, say the one on the right (it doesn't 
matter) 



2 + 



P(A)a(A) _ V'{X)^\^ 



A 



A 



1 /■« 



TT 



1 1/ 

■f A) TT A 



1 



a; 



(x — A — ze) 

(3.10) 



In the last integral use = ^ + i-K5{a) and we obtain 
P(A)c7(A) 1 



2 + 



A 



X 



(x — A) (x + A) 



(3.11) 



Let us take the derivative with respect to V{ii) (/U is > by definition of V). 

1 



a{\)5P{\) 1 d ^/\^ 



+ 



A SV{iJ.) TT 5/i /X (/Lt — A) ()U + A) 



:) 



(3.12) 



(assuming that we can show here as usual that a counterpart of {^)v,b sv{^j.) 
and (ff )v,a gy('^) vanish, see Appendix A for a proof, the exact same steps 
can be followed here). Then 



1 



(A>0) 
6p{X) 



(3.13) 



2^ vMAT^ d VW^Wo^^ 



dV{fi) 27r (7(A) dfi 
1 A/x 1 



7rV|cT(A)||a(M)|(M'-A2)^ 



A2 

[2AV - (A' +)"')(«' + ?'') + 2a V] 

(3.14) 
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Let us check immediately the b = hmit 

A 1 



(3.15) 



/RA^ Via- - A-) 
and the rest looks unfamiliar; but if we remember that we are computing 

/>2(A,//) +/32(A,-/i) (3.16) 

and 

a^-Xfi a^ + Xfi _ _ [2AV-a^(A^ + /x^)] 
(A-/i)2 + (A + /i)2 " (A2-/x2)2 ' 

(3.17) 

we check that this result agrees as expected for 6 = with the single cut 
result. Therefore for a symmetric double well, assuming no spontaneous 
symmetry breaking, we have the undisputable answer for p2{X, ^)+p2(A, —fi) 
i.e. eq. ( |3.14| ). Note that, as expected, the short distance behaviour of 
/92(A, fi) is the same as for the single well with only one cut. 

4 Asymmetric Double Well 

In order to extract p2{X,p) alone we have to consider arbitrary potentials 

instead of restricting ourselves to even (Z2) symmetric potentials as we have 
done in the above section. 
Again let us start with 

1 ^ 

pW = ]^E'^(^-^0 (4.1) 

with 

p{X) =< p(A) > . (4.2) 

Since the weight contains 

exp{-NTrV) = exp{-N'^ f ViX)p{X)dX) (4.3) 



pi(A,M) =< p(A)p(/.) >c= (4-4) 
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p(X) = -llmG{X + iO) (4.5) 
vr 

Giz) =< > (4.6) 

2G{z)=V'{z)-P{z)^) (4.7) 
aiz) = uUiz-ai). (4.8) 

See Fig. 3. 

The support of the eigenvalues consists of the two segments [ai, 02] and 
[03, 04], (we assume that they are labelled by increasing order) ; the positiv- 
ity of p{\) is satisfied provided the polynomial P{z) has an odd number of 
zeroes between 02 and 03. Contrary to the two previous cases (fl.Tj) is not 
sufficient to determine fully the polynomial P{z) and the four end points 
of the two cuts. Counting parameters and equations one sees readily that 
we miss one parameter, which we can take as the filling factor of one of 
the two wells. This factor remains undetermined at this level of the large 
N-limit, and we would have to return to a minimization of the free energy 
to fix it. However, since this parameter is not fixed at this leading order, 
we may ignore it and proceed as before for finding the leading order of the 
correlator. 



We denote 



eA = +1 as < X < 04 

-1 ai < X<a2 (4.9) 



then 

Therefore 



a{X±iO) = ±iexJ\a{X)\. (4.10) 



27rp{X)=exJ\a{X)\PiX) (4.11) 



and positivity implies that signP{X) = e\. By multiplication by 
integration over a large circle, we obtain 



(z—u){z—v) 



^ ^ P{u)a{u) - P{v)a{v) ^ V ju) ^Mu) - V (v) ^a{v) 
u — V u — V 



1 p _ r-2 dxV'{x)./\^\ 

TT Jaa Jai {x - u){x - v) ' 
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Let u = X + ie and v = fi + irj then 

^ ^ P(A)a(A)-P(/xM/x) 
A — /X 

i?(A, z/) = 

(T(A)i?(A, i^) - (7(Ai)i?(M, i^) 
A — 



<5P(A) 



-4 i^^(x)v^ZM 

as {X-X){x-Il) 



1 5 



IT du {u — X)(y — fj,) 



(4.13) 
(4.14) 
(4.15) 



(assuming once more an equivalent form of to be zero, see Appendix 

A). Hence 



a{X)R{X,iy) — = a{fi}R{n,u}- 



ir du u — X 
From these equations one finds 



(t{X)R{X,v) 



'K dv u — n 



IT du V — X 



VW)\ 



.(4.16) 



(4.17) 



and we arc left with two unknown functions /i+ and h- of a single variable. 
This gives the connected correlator 



1 



27riV2 



eAVk(A)|i?(A,M), 



I.e. 



27rArVi(A,/x) 



r^M 9 \/k(M)| ^ fee^M 1 



(4.18) 



(4.19) 



\/|o-(A)| TT M-A V|o-(//)| 

Prom its definition the two-point correlator is symmetric under exchange of 
the two eigenvalues 

P2(A,/x)=p^(/.,A). (4.20) 
This imposes the following constraints: 



h+ifj,) + h-in) = 



and 



dfj, (1 — X 



(4.21) 



(4.22) 
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A straightforward algebra gives from there the function /i+ up to an arbitrary 
constant: 



h+{X) = -{X^ -IsX + C) (4.23) 

with 

5 = 01 + 02 + 03 + 04- (4.24) 
We are thus left with one undetermined constant in the two-point function: 



47r^7V^p^(A,M) 



. ^^f^ + A^ + .^-£(A + .) + 2C). 

(4.25) 

Let us verify that, without any restriction on the constant C, this result 
satisfies the normahsation condition 



dz/p|(A, z^) = 0, (4.26) 
which fohows from the definition of P2- Returning to ( [4. 201) 



= r ^ r ^4^=0 (4^27) 

Ofi fl — X Jai Ofl fJ, — X Jaz Ofi fl — X 

since a vanishes at the end points. (This point is in fact shghtly delicate, 
since there is a non-integrable singularity at ;U = A. In the literature concern- 
ing the application of random matrices to the calculation of the fluctuations 
of the conductance in mesoscopic systems Q, this integration throughout 
the singularity is done in a routine way. A proper justification of the proce- 
dure implies returning to the true correlation function, before the smoothing 
which produces spurious short-distance singularities through replacements 
such as ( '^'"a ^ 2x^)' "^^^ smoothing is produced here by the large N 
limit.) Next consider 

(z2 - sz/2) 
dz ^ , (4.28) 

over a large circle. Note that ^17(2;) = -^^[l — ^ + 0{^)]. Consequently 
^ = \ -\- Oi\). Since there is no coefficient of - the integral vanishes. 



11 



Shrinking the contour over the cuts we obtain 

— SI// 

ex r dfi 



ir-rw-^^=o. (4.29) 



Therefore 

ra2 ra4 



ra'2 pa^ ^ no 



Again taking a very large circle 

dz 



(4.30) 

(4.31) 



since the coefficient of ^ vanishes. Therefore, shrinking the circle, 

dfi /""a dfj. 



(4.32) 



las VRmII Jai y/\cr{p)\ 

which shows that the normalization is correct for any value of C. 

Let us specialize to the symmetric double-well (ai = —b, 02 = —a, = a 
and 04 = b) 

|a(A)| = (A2-a2)(62_ a2). (4.33) 

After a few lines 



2^V2(A, M) = ,1 , ,1 , \y2 iC(^ - + A/^(V -a'- b') + a^ft^]. 

y/\a{X)\\a{p)\ [p-Xy 

(4.34) 

Several remarks (i). The result is manifestly symmetric under X <^ u 
(ii). C remains an unknown function of a & b (iii). If we assume that C 
vanishes when b=0 then since 



lim J\a{X)\ = exX{a^ - A^) (4.35) 

a — >0 * 

we recover the single band result, (iv). The cross-correlator p'^{X, —u) is sim- 
ply eq. ( 4.34 ) with u replaced by —u. Combining eq. ( [4.34[) with p^iX, —v) 
we reproduce the result eq. ( p. 14 ) in section |3[ (v). We also note that for 
C = —^{{a? + 6^) — (a + 6)^^^) we recover the connected density-density 
correlator derived from the connected Green's function (where E{k) and 



K{k) are complete elliptic integrals of first and second kind and k = f^+b) ) 
in eq. (15) of ref. p. 
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5 Orthogonal Polynomials, The Kernel, Odd and 
Even N 

This section follows the notation of ref. |^]. Let us calculate the Kernel 
KNipi, for N even and N odd 

Kn{^, = ]. (5.1) 

The asymptotic ansatz for ipn{^) for N ^ oo but N — n finite is 

V'n(A) = -^[cos{NC-{N-n^ + x + {-irv)W + 0{^)] 

/(A) = |^(^^)sin20(A) 

C'(A) = -TTpiX) 

A2 _ (a2 + ^2)/2 



(62 -a2)/2 
6 cos 4>{\) 
A 

asin</)(A) 
A ■ 

(5.2) 



cos 20(A) 

cos 2r/(A) 
sin 2r/(A) 

Substituting 

V'^(A) = ^i==[cos(7VC-(iV-iV)0 + x + (-l)^r,)(A)] (5.3) 

V^^_i(A) = ^^[cos(ArC-(iV-iV + l)</, + x + (-l)(^-^)r/)(A)] 
V /l^J 

(5.4) 

we get 

Kn{ii,v) = — ^(^==== cos iV/i(/i) cos iV/i(z^) 

X (cos cos 2ri{v) — (—1)^ sin0(i/) sin 2ri{v) 

— cos <^(/i) COS 2rj{n) + ( — 1)^ sin 4>{^) sin 2ri{ji)) 

+ sin A^/i(z^) cos Nh{fj,)(sm 4>{h') cos 2rj{i') + ( — 1)^ sin 2rj{v) cos 0(z^)) 

— sin Nh{fi) cos A^/i(i/)(sin i?^(/x) cos 27y(//) + (—1)^ sin 2r/(/Li) cos 0(/u)) 

(5.5) 
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where 

Nh{^) = iNC + x + {-l)''7j){^i). (5.6) 
On simplifying further 



[{cos{Nh{fi) + Nh{iy)) + cos{Nh{iJ,) - Nh{iy))) 

X (- cos2 - (-1)^- sin2 ^{u) - - cos2 + (-1)^- sin^ 

/i IJ, 

+ (sin(iV/i(zy) + iV/i(^)) + sin(iV/i(z>) - Nh{iJ,))) 

- {sm{Nh{n) + Nh{v)) + sm{Nh{n) - Nh{v))) 

X (- + (-1)^-) sin cos (5.7) 

Squaring and averaging we get after some tedious algebra 

Rn 



+ ^^^^(i^' + /x')(2a6(-l)^ + a2 + 62) 



2a6(-l)^aV-(a2 + 62)i^' + «'' 



2 

- {^^' + ^^')^^^^ + l{b'-ar{a' + b% (5.8) 



Simplifying for N even, 



2 RN{a + by 



X [z//x(62 + a^) - z/^/x^ - + {y - ufab] 

(5.9) 



while for N odd 



2 -RAr(a - b) 
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(5.10) 



Note that i?7V,„,„(a + hf = A{a + bf = ^^(a + bf = ^^^^ and 
R]\f^^^(a — 6)^ = B{a — b)"^ = ~ . Comparing this expression with that 
found by the previous method of section 0, we find that C = {-l)^ ab. The 
standard large N-hmit techniques of analyzing matrix models like the loop 
equation method ref. ||l], ^ and renormalization group ref. assume a 
smooth behavior with respect to N at large N. The result that C differs 
for odd or even N by terms of order one suggests that these methods may 
need to be revisited in the context of random matrix models with eigenvalue 
distributions with gaps. 

6 Conclusion 

To conclude we have outlined a method which reproduces known results 
for the single cut model and extended it to the two-cut random matrix 
model. The two-point density-density correlator contains a derivative part 
familiar from the single cut model but in addition contains a non-trivial 
non-derivative piece. It is further seen that different methods give different 
values for the two-point correlator. The orthogonal polynomial method is 
briefly outlined and gives different values for the non-derivative piece for 
even and odd eigenvalues. The loop equation method gives a different result. 
The difference in the results are in the non-derivative part of the two-point 
density-density correlator. The method outlined unifies these differences in 
a constant C which takes different values. Different values of C found from 
the orthogonal polynomial and loop equation methods are identified. 

This raises several questions regarding the analysis of this model. One 
possibility is that the even-odd differences may require some care in handling 
the large N techniques of random matrix models e.g. loop equations and 
renormalization group. Another question relates to spontaneous breaking of 
the Z2 symmetry in the large N limit. In this context, for the Z2 symmetric 
random matrix models with two wells an infinite family of solutions which 
break the Z2 symmetry and have the same free energy as the Z2 symmetric 
solution but different connected correlators have been identified in ref. [^. 
It would be interesting to compare whether the different solutions noted 
here correspond to some of the multiple solutions of ref. [^. Finally let us 
note that when the number of connected components for the support of the 
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eigenvalues changes, one finds a new universality class for the correlators. It 
is thus not completely obvious that it is legitimate to use the simple one cut 
function in the application to mesoscopic fluctuations. It seems interesting 
to us that this simple system, namely N charges confined by a symmetric 
double-well with a logarithmic repulsion between the charges, exhibits such 
rich behavior. 
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Appendix A 

Eq. (Ill) was derived under the assumption that a counterpart of 
{j§-)v,bi = 0. Here we prove this result for the asymmetric double- well. 
(Following a similar procedure we can show that an equivalent form of 
{j^)v,bi = from which eq. ( |3.12| ) follows for the symmetric double well). 
From eq. ( 4.15| ) it is easy to see that we have to prove the following equation 
equivalent to {■^)v,bi = in the single well problem 



(A-//) 



(A.l) 
(A.2) 



Ja {x - \){x - fl) IT Jc (x-A)(x-/x) 

Let us take 

2G{z) = V'iz) - P{z)^J{z - a){z - b){z - c){z - d). (A.3) 
Multiple by (^z-u)(z-v) l^z^l) ^^'^ integrate over a large circle. The first term 



Giz) dz ^^^^^ 



c {z — u){z — v) {z — a) 2m Jc z 
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as for large z G{z) ^ \, « z"^ , (z-u)(z-v) ~ 1^ and ^ « i. The 

third term becomes 

P{z)(t{z) dz 



{z — u){z — v){z — a) 27ri 

P{u){u - b){u - c){u -d)- P{v){v - b){v - c){v - d) 

{u — v) 



(A.5) 



While the second term is 



V'{z) dz _ ^^V'{u) ^ ^/^)V'{v) 



{z — u){z — v) {z — a) 2TTi {u — a){u — v) {v — a){v — u) 

^ J_ V'jx) y/^ f'i V'jx) y/^ ^ 

iri J a {x — u){x — v) (x — a) vri Jc {x — u){x — v) {x — a) 

On using u = X + ie and v = fj, + ie with A, /i on the right hand cut, 
^_x_i^ = + iTTS{x — A) and \/cr{x) = iex\y^a{x)\ the second integral 
simplifies to 

V'{z) dz 1 V'{x)\,y^ix)\dx 



{z — u){z — v) {z — a) 2-Ki tt Ja {x — X){x — fi){x — a) 

1 V'{x)\^/^\dx ^^^^^ 



TT Jc {x — A)(x — ^){x — a) 

Combining these three terms and simplifying we get eq. ( |A.2D which is what 
is needed in order to get eq. ( 4.15 ). 
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Figure captions: 



Fig. 1. The complex z-plane with one-cut and contour used for evaluat- 
ing the two-point density-density correlator for the one-cut random matrix 
model. 

Fig. 2. The complex z-planc with two-cuts and contour used for evaluat- 
ing the two-point density-density correlator for the two-cut random matrix 
model. 

Fig. 3. The complex z-plane with asymmetric two-cuts and contour used 
for evaluating the two-point density-density correlator for the asymmetric 
two-cut random matrix model. 
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Th& complex z- plane witli two asymmetiic cuts aiul coatom 
used for evaluating the two-point density-density conelator for the 
asymmetric two-cut random matrii modeL 




Fi^. 1 Tlie complex z-plane witb. one-cut and contow used for evaluAtin^ 
the two-point densitj-density conelator for the one-cut raniioni matiii modeL 




Fifi. 1 Ttss complex z-plane witli two-cuts and contour used for 
evfllufltin^ the two-point densitj-densitj coirelator fisr tbe 
two-cut random matrli model 



